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ABSTRACT 

We use a high-resolution cosmological simulation that includes hydrodynamics, multiphase star formation, 
and galactic winds to predict the distribution of metal line emission at z ~ from the intergalactic medium 
(IGM). We focus on two ultraviolet doublet transitions, O VI AA1032. 1038 and C IV AA1548, 1551. Emission 
from filaments with moderate overdensities is orders of magnitude smaller than the background, but isolated 
emission from enriched, dense regions with T ~ 10 s — 10 5 5 K and characteristic size ~ 50-100 kpc can be 
detected above the background. We show that the emission from these regions is substantially greater when 
we use the metallicities predicted by the simulation (which includes enrichment through galactic winds) than 
when we assume a uniform IGM metallicity. Luminous regions correspond to volumes that have recently been 
influenced by galactic winds. We also show that the line emission is clustered on scales ~ l/i -1 Mpc. We argue 
that although these transitions are not effective tracers of the warm-hot intergalactic medium, they do provide 
a route to study the chemical enrichment of the IGM and the physics of galactic winds. 
Subject headings: cosmology: theory - intergalactic medium - galaxies: formation - diffuse radiation 
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1. INTRODUCTION 

The modern paradigm for cosmological structure forma- 
tion, the cold dark matter model, has had admirable success 
in describing both the formation of galaxies and the distri- 
bution of matter on larger scales. The model predicts that 
small perturbations in the early universe grow through grav- 
itational instability and eventually collapse to form the ob- 
jects we see around us. As collapse proceeds, matter collects 
in a "cosmic web" of mildly overdense sheets and filaments, 
with galaxies and galaxy clusters at their intersection. Simu- 
lations and analytic studies have shown that these gas struc- 
tures are responsible for the so-called Lya forest of absorption 
lines i n quasar spectra (e . g..lBi et al.ll992tlBill993UCen et al 
1994; iZhang et alJll995t iHernauist et alJl!996t iTheuns et al 
19981 lD^vTetaiTll9 99). However, as observations and simu- 
lations have improved, it has become clear that the intergalac- 
tic medium (IGM) has experienced a variety of complex pro- 
cesses, including shock heating and feedback from galaxies. 

Most important, it now appears that a substantial frac- 
tion of the originally pristine intergalactic gas was enriched 
with a small, but significant, amount of heavy elements 
through interactions with star-forming galaxies. The best 
evidence for enrichment comes from detailed spectroscopic 
studies of the high-redshift Lya forest. Matching indi- 
vidual Lya absorption features with absorption in the cor- 
responding wavelength range for metal transitions (espe- 
cially the C IV AA1548,1551 and O VI AA1032,1038 dou- 
blets) has revealed that high-column density absorbers (with 
Nhi > 10 14 5 cm " 2 , corre sponding to physical overdensities 
of a few; e.g., Schaye 2001) ar e almost universally en- 
riched to Z > 1Q- 3 Z ? (iTvtler et alJl!995l ICowie et aUll995l 
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Songaila EUS 'Have et all EH iCarswell et al] 
2002). Pixel-by-pixel analyses have revealed the pres- 
ence of metals in even lower column density absorbers 
JCowie & Songaila! 119981 lEllison et alJ 1200(1 ISchave et alJ 
120001 120031 lAracil et alJ 120031 iPieri & Haehnelll 120031) . al- 
though the scatter in the metallicity is large and there is ev- 
idence for a decrease in m etallicity with decreasing over- 
density ( Schave et ai] l2003l) . Metals therefore appear to be 
widespread in the IGM, although the mechanism through 
which they were dispersed remains to be determined. 

The leading candidate is expulsion from star-forming galax- 
ies through winds driven b y supernovae, poss ibly aided 
by radiation pressure (e.g. Aguirre et al. 2001a[). Such 
winds hav e been observed both in nearby starburst galax- 
ies (e.g.. iLehnert & Heckmanl Il996t iHeckman et alJ UOQOt 
Martin 1999) a nd in star-forming Lyman-break g alaxies at 
z ~ 3 JPettini et al.l200lLl2002tlShaplev et alJ2003l) . but their 
power, duration, and filling factor remain unclear. There is 
also a theoretical expectation that similar winds may have 
occurred in high-redshift dwarf galaxies, because the small 
binding energies of these systems wou l d allow a wind to 
escape fairly easily dDekel & Silklll986t iMadau et"aTll2001t 
Bromm et aLll2003l [Furlanett o & Loebll2003l) . The time his- 
tory of enrichment has important implications not only for the 
characteristics of the present-day IGM but al so for such top- 
ics as the reionization history of the unive rse (IWvithe & Loebl 
I2003albt ICenll2003l ISokasian et aTll2003l) and ending the era 
of Population HI star form ation (Brom m et al.l200HlOh et ail 
l200lHMackev et alJ20 03"l). Detailed studies of the distribution 
of metals in the IGM will help to determine how and when 
enrichment occurred. Careful analysis of Lya forest absorp- 
tion spectra is one promising approach. However, absorption 
studies of the Lya forest can only provide information along 
isolated lines of sight to quasars, giving no more than an in- 
direct picture of metals in the IGM. It would thus be useful to 
find some technique to directly trace the distribution of metals 
in three dimensions. 

One such method is to search for metal line emission from 
the IGM. Because emission does not rely on the existence of 
a background source, it can be imaged over large fields of 
view. With the addition of individual line redshifts, we can 
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construct a three-dimensional picture of the metal distribution 
(in a similar way to galaxy redshift surveys). Observations 
in this direction have begun: diffuse O VI A 1032 emission 
has been detected from our own Gala xy JShelton et al.ll200H 
iDixon et alJl200U IWelsh et al]l2002b a nd above the disks of 
edge-on spirals (Ott e et alJl2Q03l) . while iHoopes et afl (2003) 
place limits on the emission from galactic winds. Because of 
the complexity of the emission processes, and their sensitive 
dependence on temperature and density, this question is best 
addressed through detailed cosmological simulations. In this 
paper, we use such a simulation to calculate emission from the 
IGM in the two ultraviolet (UV) doublets O VI AA1032, 1038 
and C IV AA1548, 155 1 . We note that doublets are particularly 
useful because they allow an unambiguous identification of 
the line feature. We focus on emission from the low-redshift 
(z < 0.25) universe that is most likely to be observable in the 
near future. We will contrast a uniform metallicity with the 
metal distribution predicted by the simulation, which includes 
winds from massive galaxies. 

Aside from enrichment, another key process affecting IGM 
gas is shock heating: simulations show that a large frac- 
tion (> 20%) of intergalactic gas has 10 s K < T < 10 7 K 
iCen & Ostrikerl 119991: iDave et alJ 12001b . This so-called 
warm-hot intergalactic medium (or WHIM) is heated through 
accretion shocks as sheets and filaments assemble into larger 
structures. WHIM gas does not form a single, well-defined 
phase: because the shocks are not uniform, it has a wide range 
of densities, from near the cosmic mean to that characteristic 
of galaxy groups. But all WHIM gas shares one crucial char- 
acteristic: because it is too warm to form stars but too cool 
to produce X-rays, it is very difficult to detect. The difficulty 
is severe enough that gas at or near the cosmic mean den- 
sity, whose surface brightness is particularly small, is referred 
to as the "missing baryons." One techniq ue to search for 
WHIM g as is to use absorption systems ilPerna & Loebi l998: 
iHellsten et alJl998b . iTripp et aiT J^OOO) found a large number 
of low-redshift O VI A 1032 absorbers along one line of sight, 
suggesting that > 10% of baryons are in the WHIM phase. 
However, these absorption studies suffer from the same draw- 
backs mentioned above and do not provide a direct picture 
of WHIM gas. Images of its distribution could thus test pre- 
dictions about the gravitational growth of structure. Simula- 
tions have shown that transitions of higher ionization stages of 
oxygen (O VII and O VIII) are e ffective probes of the h igher- 
temperature part of WHIM gas dYoshikawa et al.ll2003l) . We 
will investigate whether emission from the O VI A 1032 and 
C IV A 1548 transitions will allow one to map the distribution 
of the lower-temperature portion of the WHIM gas. Because 
gas in the appropriate temperature range is expected to he in 
sheets and filaments, this technique would offer the additional 
advantage of tracing the cosmic web, providing further con- 
firmation of the cold dark matter paradigm. 

We first calculate the line emissivity of a given gas parcel 
in ^2] We then describe our simulation and show how to cal- 
culate the emission from it in |5] We discuss our results in |4] 
and conclude in ^5] 

2. METAL EMISSION MECHANISMS 

We consider the emission from two metal line doublets: 
O VI AA1032, 103 8 and C IV AA1548, 15 5 1 . For each doublet, 
we use Cloudy 94 (Ferland 2000, 2001) to compute the emis- 
sivity of a given parcel of gas. We first construct a grid of 
densities and temperatures spanning the entire range of par- 
ticle properties in our simulation. We choose density spac- 



TableI. Ionizing Background Models 
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NOTE. — Models labeled 'HM96' and 'HM01' use the spectral shapes 
of lHaardt &M adau 1 1996) and Haardt & Madaul <2001). respectively. Model 
'Coll' assumes pure collisional ionization. All surface brightness values are 
reported in photons cm -2 s" 1 sr" 1 A" 1 . 



ings of Alog(w H /cm" 3 ) = 0.25 for 10" 6 < (n H /cm" 3 ) < 10" 2 
and Alog(«H/cm" 3 ) = 1.0 otherwise. Here «h is the to- 
tal density of neutral and ionized hydrogen. Because the 
ionization fractions (and emissivities) are sensitive to varia- 
tions in the temperature, we choose Alog(r/K) = 0.05 for 
10 4 .25 k < T < 10 7 K and Alog(J/K) = 0.25 elsewhere. We 
then generate grids of the emissivity for each transition, in- 
cluding collisional or radiative excitation followed by radia- 
tive de-excitation as well as recombination cascades. How- 
ever, photon emission following absorption of a photon of the 
same wavelength does not contribute to the net emission. We 
therefore use Cloudy to compute this "pumping" component 
separately and subtract it from the total emissivity. We gener- 
ate the grids at solar metallicity; the emissivity is proportional 
to the metallicity. 5 For each of these transitions, the shorter 
wavelength line is also the stronger member of the doublet, 
with a flux ratio of 2: 1 . For concreteness, we quote results for 
the stronger line in the following. 

Because most of the IGM is photoionized, our choice of 
ionizing background could significantly affect the results. Un- 
fortunately the ionizing backgroun d at z ~ is relatively 
poorly constrained. iDave & Trippl fcOOlb use observations 
of the low-redshift Lya forest to estimate a total H I ioniz- 
ing rate of r H i = 10" 13 - 3±0 - 7 s" 1 at z = 0.17. The shape of the 
spectrum is even more poorly known, because the extragalac- 
tic background is difficult to measure di rectly at these wave- 
lengths, particularly blueward of Lya. iBrown eta D ll2000b 
estimate a surface brightness of $(1450-1900 A) ~ 500 ± 
100 photons cm" 2 s" 1 sr" 1 A" 1 based on Hubble Space Tele- 
scope observations of diffuse radiation along a low-extinction 
line of sight through the Galaxy. The best limits in the range 
912-1216 A come from the Voyager spacecraft; Murthv et al. 
(1999) place a la upper limit of $(912-1216 A) < 
30 photons cm" 2 s" 1 sr" 1 A" 1 , although Edel stein et all ( 12000b 
argue that systematic uncertainties increase the limit to 
$(912-1216 A) < 10 4 photons cm" 2 s" 1 sr" 1 A" 1 at the 2a 
level. 

Because of this uncertainty, we include several choices 
for the radiation background. We base our models on 
the spectra presented in Haardt & Madau ( 1996, hereafter 
HM96), which includes onl y the contribution of quasars, 
and lHaardt & MadavJ ( 1200 ll hereafter HM01), which in- 
cludes both quasars and star-forming galaxies. We vary the 

5 Actually, the metallicity does affect the electron density, breaking the 
strict proportionality. However, the difference is negligible for metallicities 
in the range of interest as long as the hydrogen is highly ionized (i.e., as long 
as the gas is not self-shielded, see below). 
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FIG. 1. — The normalized emissivity of O VI A1032 . All the panels 
ionizing backgrounds (see text for details). Upper left: HM01. Upper right: 
equilibrium (solid curve). Note that this quantity is independent of density. We 
(dotted, dashed, and dot-dashed curves, respectively). 
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except the bottom right show contour plots in the n\\-T plane assuming different 
HM96. Lower left: HMOlw. Lower right: The emissivity assuming collisional 
also show s/n^ for the HM01 spectrum assuming «h = 10 -6 , 10 -4 , and 10~ 2 cm -3 



normalization s of th ese spectra over the range allowed by 
iDave & Tripd OoOl). We also include a model assuming pure 
collisional ionization and excitation. We summarize the char- 
acteristics of our ionizing backgrounds in Tabled including 
the total ionizing rate and estimates of the background levels 
at each of several wavelengths. Comparing these results to the 
direct measurements described in the previous paragraph, the 
HM01 model appears to provide a reasonable fit to the data, 
although it does underestima te the flux redward of Lya re- 
ported bv lBrown et alJ (12000). Note that, even with r H i fixed, 
the HM96 model background is much smaller than that of the 
HM01 model at the listed wavelengths. This is simply be- 
cause galaxies dominate the background redward of the Ly- 
man limit, while quasars dominate blueward of this limit. Be- 
cause it does not include galaxies, the HM96 model strongly 
underestimates the background at ~ 1000-2000 A. 

The resulting emissivity grids are shown in Figure ^ f° r 
O VI A1032 and Figure 0for C IV A1548. We display, using 
logarithmic contours, the density-normalized line emissivity 
ejn\ (in units of erg cm 3 s" 1 ). In such units, the emissivity 
in a purely collisionally ionized gas is independent of density. 
Three of the panels show contour plots of the emissivity for 
different choices of the ionizing spectrum: HM01 (top left), 
HM96 (top right), and HMOlw (bottom left). We draw con- 
tours at intervals of Alog(e/«^) = 1.0 and label the highest 
contour in each panel for clarity. The bottom right panel in 
each figure shows a line plot of e/n^ as a function of tem- 
perature. The solid line assumes pure collisional ionization, 



while the dotted, dashed and dot-dashed curves assume the 
HM01 background with n H = 10" 6 , 10" 4 , and 10" 2 cm -3 , re- 
spectively. We see that O VI A 1032 efficiently probes gas with 
T ~ 10 5 5 K, while C IV A1548 probes gas with T ~ 10 5 K. 
Note that the line emissivities are quite sensitive to the tem- 
perature, and a realistic prediction of the IGM emission re- 
quires an accurate treatment of heating and cooling processes. 
The bottom right panels show that photoionization has a sub- 
stantial effect on the emissivity of the low-density IGM: it 
increases the emission from cool gas but decreases the emis- 
sion from hot gas. The photoionizing background drives low- 
density gas to higher ionization states than it would other- 
wise assume. On the other hand, these panels also show 
that the photoionizing flux makes little difference for gas with 
«h/«h ^ 100. Such gas is dense enough for collisional pro- 
cesses to dominate over any reasonable ionizing background. 
A comparison of the HM01 and HMOlw panels shows that, 
for lower densities, the emissivity quickly returns to the colli- 
sional ionization values as the ionizing photon flux decreases 
(because the number of available ionizing photons decreases, 
so the maximum gas density that can be affected declines). 
Finally, note that the HM96 and HM01 spectra yield nearly 
identical results, despite the large differences in the back- 
ground fluxes for A ~ 1000-2000 A. This is because we have 
removed the pumping contribution to the emissivity and be- 
cause the ionization potentials for these species are high, well 
into the regime in which quasars dominate the spectrum. 
One uncertainty in our emissivity calculations is that they 



FURLANETTO ET AL. 




4.5 5.0 5.5 6.0 6.5 4.5 5.0 5.5 6.0 6.5 

log(T/K) log(T/K) 




4.5 5.0 5.5 6.0 6.5 4.5 5.0 5.5 6.0 6.5 

log(T/K) log(T/K) 



FIG. 2.— Same as FigureQ] but for C IV A1548. 



implicitly assume ionization equilibrium. Testing this as- 
sumption precisely is difficult because equilibrium for high- 
ionization state metals depends on a large number of separate 
recombinations. However, it appears to be a reasonable as- 
sumption. We have used Cloudy to compute the recombina- 
tion (and collisional ionization) times of C IV and O VI, along 
with those of the neighboring ionization states, in the range 
4.75 < log(r/K) < 6.0. For metallicities Z < 0.1Z Q , the re- 
combination and ionization timescales are generally smaller 
than the cooling timescale by a factor of a few. For signifi- 
cantly higher metallicities, the times become comparable. In 
this case, the emissivity of gas near the excitation tempera- 
ture of a given transition will decrease (because the relevant 
ionization state would be less abundant than expected in equi- 
librium). However, the gas will emit over a wider temperature 
regime - and hence for a longer time - as it passes through the 
O VI and C IV states. The net effect will be to make the emis- 
sivity somewhat less sensitive to the gas temperature, but only 
in the highest metallicity regions. 

We also consider two hydrogen lines that may be sources of 
background contamination in metal line studies: Lya A1216 
(which may be confused with C IV A 1548 emitted from lower 
redshifts) and Ly/3 A 1026 (which may be confused with O VI 
A1032). If an instrument can resolve the O VI or C IV dou- 
blets, these lines may be unambiguously identified. Other- 
wise, single spectral lines like the Lyman series will mimic 
metal emission. We calculate the total emission in both Lya 
and Ly/3 using Cloudy, with the same grid spacings as above. 
Subtracting the pumping contribution is complicated by the 
fact that excitations from absorbed background photons do 



not necessarily produce a Ly/3 photon: the excited state can 
instead decay to n = 2 and then to the ground state. (For ex- 
ample, a Ly/3 photon can be "converted" into Lya and Ha 
photons.) In this case the number of photons removed from 
the background radiation field by Ly/3 does not necessarily 
equal the number emitted. We compute the amount of ab- 
sorption through 

Epump = 2im m hvB n J, (1) 

where «hi is the neutral hydrogen density, B\2 is the Einstein 
absorption coefficient, v is the frequency of the line, and / 
is the background flux averaged over the line profile. For a 
reasonably smooth background spectrum, J is simply equal 
to the background flux evaluated at the line center. We sub- 
tract the amount of energy absorbed from the total emissivity 
to compute the net emissivity. Note that the net result can be 
negative if photon conversion dominates over collisional ex- 
citation and recombination; in this case we set the emissivity 
to zero. We compute and subtract the pumping contribution to 
Lya using C loudy, just as for the m etal transitions. We refer 
the reader to Furlan etto et alJ J2003I) for more information on 
Lya emission. 

3. THE COSMOLOGICAL SIMULATION 

3.1. Simulation Parameters 

We perform all of our calculations within the G5 cosmo- 
logical simulation of ISoringel & Hernauistl fe003bl) . This 
smoothed-particle hydrodynamics (SPH) simulation uses 
a mod ified version of the GADGET code ( Spring e! et alJ 
2001b) incorporating a new conservative formulation of 
SPH with the specific entropy as an independent variable 
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llSpringel & HernauisJ 120021) . This "conservative entropy" 
approach has several advantages over conventional treatments 
of SPH. Since the energy equation is written with the en- 
tropy as the independent thermodynamic variable, the ' pdV 
term is not evaluated explicitly, reducing noise. Including 
terms involving derivatives of the density with respect to the 
particle smoothing lengths enables this approach to conserve 
entropy (in regions without shocks), even when smoothing 
leng ths evolve adaptively, avoiding the problems noted by, 
e.g., [Hernauist (1993), and moderates the overcooling prob- 
lem present in earlier formulations of SPH. 

The simulation assumes a ACDM cosmology with fi m = 
0.3, n A = 0.7, n h = 0.04, H = 100/i km s" 1 Mpc" 1 (with 
h = 0.7), and a scale-invariant primordial power spectrum with 
index n = 1 normalized to a% = 0.9 at the present day. These 
parameters are consistent with the mo st recent cosmological 
observations (e.g., Sperg eTet al.ll2003T) . The particular simu- 
lation we choose has 324 3 dark matter particles and (initially) 
324 3 SPH particles in a box with sides of 100/T 1 comov- 
ing Mpc, yielding particle masses of 2.12 x 10 9 /z _1 M Q and 
3.26 x 10 h~ l M Q for the dark matter and gas components, 
respectively. The spatial resolution is 8/T 1 kpc (comoving). 

Figures [0 and |2] demonstrate that the emission character- 
istics of the IGM are extremely sensitive to its temperature 
and density. We show in Figure [5] the distribution of simu- 
lation particles (at z = 0) in the n^-T plane. Note that the 
colorscale is logarithmic. The phase diagram shows several 
distinct loci of points. The majority of particles lies along a 
narrow line connecting underdense, cool gas to moderately 
overdense, warm gas. This curve represents the unshocked 
IGM gas, and its shape is fixe d by a c ombination of photoheat- 
ing an d adiabati c expans ion ( Katz et al.ll 1 9961 iHu i & Gnedin 
ll997USchave et alll999HMcDonald et all2001l) . The second 
set of particles occupies a much broader range in tempera- 
tures at moderate to high overdensities. It represents shocked 
gas in filaments or virialized halos, including the WHIM gas. 
A third set of particles lies along a nearly vertical line at 
T ~ 10 4 K, representing gas that has cooled and collapsed 
into bound objects. The rapid decrease in the cooling rate at 
T < 10 4 K forces such gas to approach a constant asymptotic 
temperature slightly below that value. Finally, a fourth set of 
particles lies at high densities («h ^0.1 cm" 3 ) and includes 
all particles able to form stars. 

The simulation includes a new method to describe star 
formation a nd feedback in the interstell ar medium (ISM) 
of galaxies ( Springel & Hernquist 2003a). The model di- 
vides the ISM into a cold, star-forming component and a 
low density component heated by supernovae (stars are also 
included separately). Gas can move between the phases 
through three processes: star formation, the evaporation of 
clouds by supernovae, and cooling of the hot phase. This 
model yields a clearly defined density threshold for star for- 
mation («h > 0.129 cm" 3 ). Most important, it yields a con- 
verged numerical prediction for the cosmic star formation 
rate within the extended serie s of simulations performed by 
Snrin gd"& Hernau isj ll200 3b|). In fact, usi ng simple physi- 
cal arguments, Hernauist & Springel (2003) have shown that 
the evolution of the star formation rate density in these sim- 
ulations has a simple analytic form that is determined by the 
gravitational growth of halos (the limiting factor at high red- 
shifts, when cooling times are short) and the expansion rate 
of the universe (the limiting factor at low redshifts, when 
cooling times are long). W e refer the interested reader to 
lHernquist"& Springel (2003) for more details on the mech- 
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FIG. 3. — Phase diagram of simulation particles in the box at z = 0. The 
colorscale represents the relative particle density (in logarithmic units). The 
horizontal solid line marks the threshold density for star formation; we ex- 
clude all particles above this line from our analysis. Similarly, the dashed 
box marks particles expected to be associated with self-shielded regions (also 
excluded from our analysis). 



anisms driving the star formation histor y. The sui te of simu- 
lations performed by Springel & Hernauist (2003b) together 
with the converged star formation prediction give us confi- 
dence that total star formation rates in the simulation we have 
chosen are unaffected by resolution during the era we study 
(z<0.5). 

The star formation model determines the location of star- 
forming particles in Figure [3] We show the threshold den- 
sity for star formation by a horizontal solid line. For parti- 
cles above this threshold, the simulation reports the average 
temperature of the cold and hot phases. Dense, actively star- 
forming gas populates the narrow locus of particles with mean 
temperatures reaching ~ 10 6 K, while the faint vertical strip 
at 10 4 K marks particles that have just been selected by the 
phenomenological wind model to be ejected from the star- 
forming ph ase (see below). Using the method described in 
Sprin gel & Hernauist! l l2003al) . we could in principle separate 
the cold (star-forming) gas and the hot (supernova-heated) gas 
in order to predict the metal line emission from each phase. 
This would be appropriate if we wished to compute the emis- 
sion from galaxies themselves. However, such a treatment 
would ignore the inevitable uncertainties in the local ioniz- 
ing radiation field (which would be much greater than our as- 
sumed diffuse background), dust content, and the geometry of 
the two phases within each galaxy. Given that the simulation 
ignores all of these processes, we chose to simply eliminate all 
star-forming particles from our analysis. As a consequence, 
our results should be interpreted as excluding emission from 
galaxies. Galaxies will actually emit strongly in the lines we 
study (at least relative to pixels without galaxies), but they 
provide little information about the IGM and can be ignored 
for the purposes of this investigation. 

We also exclude self-shielded gas from our analysis: by 
definition, the ionizing radiation field in such regions is much 
smaller than the mean field. We would therefore overesti- 
mate the ef fects of photoion ization in such regions . Both 
simul ations dKatz et all 19 96) and analytic estimates (Schave 
2001) show that self-shielding will only become important 
in dense, cool regions. We therefore exclude all gas with 
«h > 10~ 3 cm" 3 and T < 10 4 5 K from our analysis. This den- 
sity threshold marks the point at which the neutral hydrogen 
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column density Nm ~ 10 17 cm -2 (e.g., Sch avdl2001l) . corre- 
sponding to an optical depth to ionizing radiation of unity. 
These criteria are shown by the dashed lines in Figure |3j we 
see that self-shielding is only important for gas on or near 
the cooling locus described above. A comparison with Fig- 
ures \l\ an d Q show that excluding self-shielded gas has only 
a small effect for C IV A1548 and O VI A1032: the emissiv- 
ity of gas below this temperature threshold is so small that 
our results are robust to uncertainties in self-shielding. How- 
ever, gas in this region of phase space will produce copious 
numbers of Lya photons. Our results in this paper there- 
fore underestimate the number of bright Lya emitting regions; 
iFurlanetto et alJ (120031) consider Lya emission in much more 
detail. 

Note that we do include particles that do not actively form 
stars but nevertheless lie near galaxies. Of course, local 
sources of ionizing radiation will also affect such particles. 
However, bec ause bright particles are primarily collisionally 
ionized (see §4.6> . we do not consider this approximation a 
bad one. Figure^shows that photoionization becomes impor- 
tant when the local ionizing radiation field exceeds the diffuse 
background by a factor ~ 100. This condition will only be 
satisfied near galaxies; including local ionizing sources will 
probably increase the signal because much of the gas near 
galaxies is cooler than the excitation temperature of these 
transitions. A more serious concern is dust: because many 
of these particles have recently been ejected from galaxies, 
they may have significant amounts of dust that redistributes 
flux emitted through metal lines. We do not attempt to model 
this process because of the uncertainties involved. 

The simulation also includes a model for galactic winds. 
Such winds are the most plausible mechanism for carrying 
metals from the star-forming regions in which they are cre- 
ated to the low-density IGM i n which they have now been se - 
curely detected (e.g.. lElrison et alJl2000tlSchave et alll2003l) . 
The volume over which these winds penetrate is clearly crit- 
ical for understanding the distribution of metals (and hence 
the distribution of O VI A 1032 and C IV A 1548 emission) 
so we review the wind model here (see Springel & Hernquist 
20033 for more details). Because of the complexity of the 
wind generation mechanism in individual galaxies, the treat- 
ment must necessarily be phenomenological. We first assume 
that the disk mass-loss rate in the wind, M w , is proportional 
to the star formation rate : M w = r)M±. Following the obser- 
vations of Martin ( 1999), we set 77 = 2. We further assume 
that essentially all of the available supernova energy powers 
the wind. These conditions fix the wind's initial velocity at 
v w = 484 km s" 1 . Each particle in a star-forming region is as- 
signed a probability based on M w of escaping from the galaxy, 
and those particles that escape are given a velocity boost in a 
random direction in such a way that momentum is statistically 
conserved. Particles that have recently been selected for ejec- 
tion appear as a high-density "tail" to the cooling IGM gas in 
Figure y] It is worth noting that we exclude emission from 
these particles, although in reality they will contain some en- 
riched, warm gas that may produce metal line emission. How- 
ever, these particles are directly associated with galaxies and 
would be difficult to distinguish from their host galaxies. 

Note that the simulation itself included a diffuse ionizing 
backg round with the shape of HM96, normalized as described 
in lDave et alJ (119991) to match Lya forest observations. This 
choice leads to reionization at z — 6 and (after this time) 
largely determines the temperature distribution of the diffuse 



gas. We therefore cannot include the different heating effects 
of the ionizing backgrounds described in Table 1 ; fortunately 
this will have only a small effect on the observable results. 
A much more important uncertainty is the cooling function 
used in the simulations, which does not include metal line 
emission. To estimate the importance of this effect, we used 
Cloudy to compare the cooling times (including metal lines) 
of gas in the range 5 < log(T/K) < 6 for several different 
metallicities. We find that for Z = O.1Z , the metal lines do 
have a substantial effect on the cooling times, decreasing them 
by a factor of 2-4 over metal-free gas. However, the cooling 
times fall by only 10-30% for Z = 0.01Z©. Because the latter 
value is closer to the mean metallicity in the simulation, we 
do not expect metal line cooling to have a significant effect 
in most of the IGM. However, it will affect the characteris- 
tics of highly enriched gas, such as that in winds. We would 
therefore expect winds to cool more quickly through the tem- 
perature regime in which they emit strongly, decreasing the 
number density of such sources by a factor of a few. At the 
same time, each individual source would be somewhat more 
compact and hence brighter when it passed through the rel- 
evant temperature range. However, we expect the real effect 
of metal line cooling to be somewhat smaller than this esti- 
mate for two reasons. First, increasing the cooling rates will 
also increase the amount of gas able to form stars, which will 
in turn increase the amount of energy put into winds. Sec- 
ond, the above cooling time estimate does not include cooling 
due to the expansion of t he wind. In semi-analytic mod els of 
winds (e.g. Ma dau et all200 It IFurlanetto & Loebl2003l) . this 
cooling mechanism is quite important, so the actual differ- 
ence between the cooling times will be smaller than suggested 
above. 

3.2. Processing the Simulation 

Given the emissivity grids and the simulation outputs, we 
now describe how to calculate the metal line emission from 
a slice of the universe at redshift z with a specified depth Az 
and angular width 9 S . We first identify the simulation output 
(or outputs) containing this slice. (Note that the outputs are 
spaced one light-crossing time apart, so interpolation between 
them is negligible; see Snringel et al. 2001a). In order to en- 
sure that we sample a random volume of the box, we perform 
the fo llowing steps before processing the output dCroft et all 
120011) . For each box, we apply a random translation along 
each axis, randomly choose an axis of the simulation box to 
be the line of sight, and randomly reflect around each axis. 
If 6 S < Obox/ V2, we also rotate each particle around the line 
of sight by a randomly generated angle (when this criterion is 
not satisfied such a rotation would interfere with the periodic 
boundary conditions imposed on the box). We then divide the 
projected volume into a grid of A^ 2 pixels. In the results shown 
here, we choose N = 600. 

Once these steps are complete, we compute the emission on 
a particle-by-particle basis. For each particle in the box, we 
first check whether it lies within the slice of interest. For this 
purpose, we include peculiar velocities in assigning the parti- 
cle a redshift. We then assign to each particle within the slice 
a volume V; =M,/p/, where M; is its mass and p,- is its mean 
density. For simplicity, we assume that each particle is cubical 
with a uniform density. In principle, we should distribute the 
particle mass and density following the SPH kernel. However, 
we find that luminous particles are typically smaller than our 
map resolution, so smoothing would make little difference to 
our conclusions. Low-density particles in the diffuse IGM are 
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FIG. 4. — Maps of a 2.5° X 2.5° slice of the universe with z = 0.07, Az = 0.01, and AO = 1'. The top row shows the surface brightness of the C IV A1548 
(upper left) and O VI A1032 (upper right) transitions assuming simulation metallicities. The colorscale is labeled in units of photons cirT 2 s _I sr -1 . The bottom 
row shows the projected mass-averaged temperature (lower left) and metallicity (lower right). 



often larger than our grid cells, but their emissivity is signif- 
icantly below any realistic detection thresholds (see below). 
We have checked that the particle widths given by our proce- 
dure are roughly consistent with the SPH smoothing lengths 
determined from the simulation. 

We then assign a metallicity to each particle. We use two 
prescriptions that contrast two different enrichment processes. 
One option is to use the metallicities reported by the simu- 
lation. These are ultimatel y de termined by the wind feed- 
back process described in 33.11 and thus reflect the enrich- 
ment pattern expected of powerful winds from massive galax- 
ies. Alternatively, we can assume that metals are distributed 
uniformly throughout the IGM. In this case we normalize the 
total mass of metals throughout the output box to that of the 
simulation. We note that, in our simulation, the mean mass- 
weighted metallicity of the gas at z = is Z = 0.024Z Q . Of 
course, for a uniform distribution the emissivity of each pixel 
is proportional to the mean metallicity of the simulation and 
can be res eated to any desired level. We use the solar abun- 
dances of lAnders & Grevessel dl989l) : (C/H) Q = -3.45 and 
(O/H) = -3.13. 



We determine the emissivity of each particle from its den- 
sity and temperature by linearly interpolating (in log space) 
the appropriate grid described in ^2] and multiplying by the 
metallicity of the particle (in solar units). The surface bright- 
ness $, of the z'th particle (in units of photons cm" 2 s" 1 sr" 1 ) 
is then 

$ . = !+£■._£*_ JL 
' 47rr2 ,. 27rfa/oAfi;' 

where vq is the intrinsic frequency of the line, ru is the lu- 
minosity distance to the particle, e, its emissivity, Zi its red- 
shift, and AO; the solid angle subtended by the particle. We 
distribute this flux over pixels on the sky by computing the 
fraction of the particle volume that overlaps each pixel. The 
total surface brightness of a pixel j is then 

$C/) = ^$. A « (3) 

i 

where Aj is the fraction of pixel j subtended by particle i and 
the sum is over all particles. Finally, we smooth the maps us- 
ing a Gaussian filter with a width of four pixels; the resolution 
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A9 quoted below is the FWHM of this smoothing function. 



4. RESULTS 
4.1. Metal Line Maps: Qualitative Features 

Figure El shows maps of the C IV A1548 and O VI A1032 
surface brightness (upper left and upper right, respectively, 
both assuming simulation metallicities) together with maps 
of the projected mass-averaged temperature (lower left) and 
metallicity (lower right) within a slice at z = 0.07. The red- 
shift width is Az = 0.01, corresponding to AA ~ 10 Afor O VI 
A1032; a high spectral resolution instrument would therefore 
be able to probe considerably narrower redshift slices. 6 All 
the maps have AO = 1' (« 55/i _1 physical kpc at this redshift) 
and use the metallicities reported by the simulation. Sev- 
eral points are immediately apparent. First, the two transi- 
tions trace roughly the same structures: it is rare for a re- 
gion to emit in one without the other (although this is not 
necessarily true on a pixel-by-pixel basis), because they both 
trace warm high-density gas (with T ~ 10 K for C IV and 
T ~ 10 5 5 K for O VI). Filaments are warmer, denser, and 
more enriched than the surrounding voids so they stand out 
clearly in both transitions. This is a simple consequence of 
the metal dispersal mechanism: winds must begin in galaxies 
which in turn lie inside overdense regions. (Note, however, 
that the metallicity is quite inhomogeneous even within fila- 
ments.) Second, the O VI maps better trace the diffuse fila- 
mentary structure of the IGM. In the C IV maps, filaments ap- 
pear to be composed of discrete clumps rather than a smooth 
gas distribution; this transition traces cooler gas that is pre- 
sumably more tightly associated with galaxies. Finally - and 
most important - the mean surface brightness of these tran- 
sitions is extremely small, with O VI A 1032 is significantly 
stronger than C IV A 1548 . Given the expected diffuse back- 
ground (> 10 photons cm" 2 s" 1 sr" 1 A" 1 for O VI A1032 and 
> 100 photons cm -2 s" 1 sr -1 A -1 for C IV A1548; see Table 
1), only the brightest regions will be visible. 

Figure|4]shows that the brightest metal-emitting regions are 
associated with hot, approximately spherical regions strung 
along filaments. These regions are heated by two mecha- 
nisms. One is structure formation: T ~ 10 6 K gas is charac- 
teristic of large galaxies and small galactic groups. The large 
hotspot in the upper left of this map is such a group. A sec- 
ond heat source is galactic winds. In our feedback model, 
the characteristic temperature of the wind is expected to be 
Twind ~5x 10 6 /heat K, where /heat is the fraction of the ini- 
tial kinetic energy of the wind that remains in thermal energy. 
Of course, /h ea t is difficult to calibrate because it depends on 
cooling within the wind medium (including adiabatic cool- 
ing as the wind expands) and on the gravitational potential 
from which the wind has to escape, but ongoing winds can 
clearly heat gas to the range relevant for these two transitions. 
Winds like the ones in our prescriptio n can reach distance s 
~ 1/T 1 Mpc from their host galaxies ( Agui rre"et alJ l2001b). 
the scale of many of the smaller, nearly spherical hotspots in 
Figure^ 

We can better understand the qualitative features of these 
maps by looking more closely at which particles are lumi- 

6 Note that we ha ve chosen slice s that are 10 times thicker (in redshift 
space) than those of Furlanetto et al. 1 2003), because the number density of 
metal emitters is much smaller than that of Lyct emitters. 



nous. Figure |5] again shows the distribution of particles in 
the rift-T plane, but this time we have weighted the number 
density by the emission in the separate transitions. The top 
row weights the density by O VI A 1032 emission assuming 
uniform (left panel) and simulated (right panel) metallicities. 
The bottom row weights the density by C IV A 1548 emissiv- 
ity. Note the considerable stretch in the colorscale; gas in 
all but the darkest regions will in fact be unobservable (espe- 
cially considering that the low-density gas is distributed over 
a much larger volume). We see in the figure that, in all cases, 
the total emission comes primarily from a small region within 
the plane: dense gas in the range T ~ 10 s — 10 6 K. This is re- 
assuring in two ways. First, a comparison to Figures ^ an d 
[2] shows that the gas responsible for most of the emission is 
dense enough for its emissivity to be essentially unaffected 
by the choice of photoionizing background. We thus expect 
our results to be robust against the uncertainties in the ion- 
izing background. Second, note that O VI A 1032 emission 
from dense gas vanishes for T < 10 5 K. Because this is well 
above our self-shielding threshold (see Figure[3), the location 
of this cut will have little effect on our predictions about this 
transition. C IV A 1548 emission does extend below the self- 
shielding threshold (the cut is most visible through the hori- 
zontal line at «h = 10~ 3 cm -3 ). However, most self-shielded 
gas rapidly approaches T ~ 10 4 K as it becomes denser; such 
temperatures are well below the ionization threshold of C IV 
A 1548 and hence the actual emission from self-shielded gas 
turns out to be small. Furthermore, note that the emission 
from cooling gas is quite weak in the model with simula- 
tion metallicities, so we would not expect self-shielded gas 
to make as much of a difference in this case. 

However, the news is not all positive: a comparison to Fig- 
ure [3] shows that most of the emission comes from gas lying 
between the cooling and IGM loci, especially in the simula- 
tion metallicity model. The area in which gas emits strongly 
is only weakly populated in the simulation and hence metal 
emission is concentrated within rare, bright regions. We also 
conclude that these transitions are not efficient probes of gas 
on the standard IGM n^-T curve or of the shocked WHIM 
gas. Finally, comparison of the two different metallicity mod- 
els shows that a larger fraction of emission comes from gas on 
the IGM and cooling loci if we assume a uniform metallicity. 
This is a simple consequence of the metal distribution illus- 
trated in Figure|4] metals tend to be more concentrated inside 
filaments, which are already relatively bright emitters. A uni- 
form metallicity thus decreases the contrast between filaments 
and voids. 

As mentioned above, nearly all of the emission in Figure 
[4] is much less than the background and hence unobserv- 
able. Figure [6] shows an example of what one could see in 
a mock (background-limited) observation. Panels (a) through 
(e) show a slice of the universe with z = 0.15, Az = 0.01, and 
AO = 0.24' (corresponding to a physical size « 25/T 1 kpc). 
Figures I6J3 and [(J> show the surface brightness of the O VI 
A 1032 transition assuming uniform and simulation metallici- 
ties, respectively. Figures |6j/ and |6^ show the same for the 
C IV A 1548 transition. We include only those pixels with 
$ > 7.5 photons cm" 2 s -1 sr" 1 . Table 1 shows that this is a 
reasonable guess for the diffuse background for O VI A 1032, 
assuming a spectral re solution < 1 A. It i s also compatible 
with the upper limit of Murth v et alJ {l299). However, it un- 
derestimates, by at least an order of mag nitude, the back- 
ground redward of Lya (Brown et al. 2000). It will therefore 
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FIG. 5. — Distribution of simulation particles in the box at z = 0. Each panel shows the relative number density of particles in the na-T plane, with a particular 
weighting. Note that the contours are logarithmic. The top row weights the number density by O VI A1032 emissivity for uniform (upper left) and simulation 
(upper right) metallicities. The bottom row does the same for the C IV A1548 transition. We have excluded star-forming and self-shielded particles. 



be extremely difficult to pick out the weaker features shown in 
the C IV A1548 maps. In all cases, we see that potentially ob- 
servable emission is confined to isolated bright regions. The 
number density of such regions is small, particularly for C IV 
A1548, even though we have chosen a rich region for illus- 
trative purposes. This slice has about three times more bright 
O VI A 1032 regions and six t imes more bright C IV A 1548 
regions than average (see §4.3> . As expected, the peak bright- 
ness increases if we use the simulation metallicities, because 
these regions are much more enriched than average. We find 
that the typical sizes of the bright spots are ~ 80-150/T 1 phys- 
ical kpc for O VI A1032 and - 40-50/T 1 kpc for C IV A1548. 
Clearly, O VI A 1032 emission reveals much more interesting 
structure than C IV A 1548 can. 

Figure shows another example of a mock observation. 
The left panel shows a 1° x 1° map of the O VI A 1032 transi- 
tion with AO = 0.24' and simulation metallicities. We show 
a full "light cone" from z = 0.001 to z = 0.16, illustrating 
what one could observe with a wide-field integral field spec- 
trometer, and we again include only those pixels with <£> > 
7.5 photons cm" 2 s" 1 sr" 1 . The right panel shows a closeup of 
a rich section of the larger map. A typical pointing of such an 
instrument would detect a large number of emission regions. 
It is also obvious that bright regions are rarely spherical: many 
features have complex substructure that reflects their internal 
dynamics. 

4.2. Contaminants 

As discussed in |2j emission from neutral hydrogen is a 
potential contaminant in searches for these lines. Figure [(J- 



shows the Ly/3 emission from the same slice as the other 
panels while Figure^ shows a different volume correspond- 
ing to the same observed wavelength interval as the C IV 
A1548 panels (placing it at z = 0.47). 7 Emission from the 
Ly/3 line is morphologically similar to that from C IV A1548, 
though it is even more concentrated. Strong Lya or Ly/3 
emission requires a relatively neutral, dense medium and so 
comes primarily from gas that is near to being self-shielded. 
The total amount of emission in these lines would have been 
markedly greater if we had included self-shielded gas, but 
it would have re mained concentrated in isolated objects (see 
Furlanetto et al. 2003, which contains a detailed discussion of 
Lya emission). The probability of a chance superposition of 
bright C IV A1548 and Lya emission remains small. 

4.3. Pixel Probability Distribution Functions 

Figure [8] shows surface brightness information about our 
maps in a more quantitative form. We have generated an 
ensemble of 100 maps, each containing 150 2 pixels with 
AO = 0.24', and computed the probability distribution func- 
tion (PDF) of pixel surface brightness. Here the PDF is de- 
fined as dn/dln^ (with unit normalization). The top and bot- 
tom panels show the PDFs of the O VI A 1 032 and C IV A 1 548 
transitions, respectively. In each panel, we show results 

7 Actually, for a given observed wavelength, Ly/3 and O VI A 1032 also 
trace different volumes. However, the wavelength difference between the 
two transitions is sufficiently small that we show the same volume for ease 
of comparison in the emission mechanisms. The statistical properties of Ly/3 
emission are unaffected by the small redshift translation. 
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assuming both simulated (solid lines) and uniform (dashed z = 0.47 (bottom panel), illustrating the expected background 
lines) metallicities. The dotted curves show the PDF for the from these transitions. The insets show the bright end of the 
surface brightness of Ly/3 at z = 0.15 (top panel) and Lya at pixel distribution. We remind the reader that these histograms 
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Ly«, see text), Az = 0.01, and A9 = 0.24'. The upper and lower panels show 
the PDFs for the O VI A1032 and C IV A1548 transitions. In each panel, the 
solid and dashed curves correspond to simulation and uniform metallicities, 
respectively. The dotted curves show the background from Ly/3 and Lya 
transitions (upper and lower panels, respectively). The insets show the bright 
end of the PDF. All curves assume the HM01 ionizing background. 



show only the O VI A 1032 and C IV A 1548 components of the 
doublets; the total emission from the doublet is 50% greater 
than the numbers we quote. On the other hand, if both mem- 
bers of the doublet are required to confirm the reality of a 
feature, the weaker member fixes the experiment's sensitivity. 

Bright regions are obviously rare (particularly for C IV 
A 1548), covering only about one percent of the sky in 
this redshift range. Because the sources are extended, the 
PDFs do not translate directly into space densities. We 
find that the comoving density of distinct sources with $ > 
100 photons cm" 2 s" 1 sr -1 is «ovi ~ 7 x 10 _3 /i 3 Mpc -3 and 
"civ ~ 3 x 10 _3 /z 3 Mpc~ 3 (assuming simulation metallicities); 
the difference between the two is less dramatic than suggested 
by the PDF because O VI A 1032 regions are much more ex- 
tended. Clearly, metallicities from the simulation decrease 
the surface brightness of the low-density IGM but dramati- 
cally increase the number of bright pixels. Luminous regions 
have recently been enriched by galactic winds and hence have 
above average metallicity; on the other hand, simulation parti- 
cles at or below the mean density of the universe have below- 
average metallicities (because they are far from the galax- 
ies that enrich the IGM). We also see that PDFs calculated 
from simulation metallicities have broad tails extending to- 
ward high surface brightness; these correspond to enriched 
regions that are moderately overdense and/or outside of the 
peak temperature range of the transition (i.e., filamentary gas 
in Figure |6). 

Moreover, we find that Ly/3 emission from the low-density 
IGM is negligible, even compared to that from O VI A 1032, 
but that Lya emission from this gas dominates the metal lines. 
The difference between the two hydrogen lines occurs be- 
cause (1) a much larger fraction of recombinations result in 
Lya photons than Ly/3 and (2) the absorption of a A = 1026 A 
photon by a neutral hydrogen atom need not result in re- 
emission of a Ly/3 photon; it can instead be converted to Ha 
and Lya photons. In the diffuse IGM, this conversion process 



dominates over the Ly/3 photons produced through recombi- 
nations and collisional excitation, so the net emissivity is neg- 
ligible. The inset shows that pixels are more likely to host 
Lya emission than C IV A 1548 emission, so having the spec- 
tral resolution to resolve both members of each doublet will be 
crucial to UV metal line emission studies. Contamination by 
Lya emission becomes eve n more severe if we inc lude emis- 
sion from self-shielded gas dFurlanetto et al 12001 . 

As a check on our methods, we now estimate the mean sur- 
face brightness of the Lya transition. We first note that (ex- 
cept in self-shielded regions) the gas is highly photoionized; 
in this case recombinations dominate over other Lya pro- 
duction processes. The recombination rate is n rec = cmn e n p , 
where an is the recombination coefficient, and a fraction 
r\ « 0.42 of recombinations produce Lya photons (~ 38% of 
recombinations go directly to the ground state, while ~ 32% 
of recombination s to excited states p ass through the 2 S rather 
than the 2 P state; Osterbrock 198 9j). The surface bri ghtness 
per unit wavelength is then (Gould & Weinberg 1996) 
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where h(z) = H(z)/Hq and 6 represents the effective overden- 
sity of the observed pixel. Given the bandwidth of ~ 15 A 
in Figure |8j this simple argument reproduces the peak of the 
Lya PDF with S ~ 1, giving us confidence in our analysis 
procedure. 8 We note that this formula only applies to opti- 
cally thin absorbing clouds and cannot describe bright pix- 
els in or near the self-shielded regime JGould & Weinberg l 
1996). In principle, a generalization of this formula could be 
used to estimate the metal line signals as well. However, the 
temperature-dependent ionization states of oxygen and carbon 
significantly complicate any such attempt; we can really only 
extract from it a rough guide to the expected redshift evolution 
of metal line emission from the low-density IGM. 

4.4. Angular Resolution 

As mentioned above, O VI A 1032 regions have physical 
scales ~ 80-150/i _1 kpc and C IV A1548 regions have sizes 
< 40/T 1 kpc. Figures [6] and show that emitting regions 
can have complex substructure. Typically, O VI A1032 emit- 
ters have bright ($ > 500 photons cm -2 s _1 sr -1 ) central re- 
gions of size > 40/T 1 kpc surrounded by irregular, low sur- 
face brightness halos. C IV A1548 regions lack the surround- 
ing halos. Resolving this substructure requires reasonably 
good angular resolution. We show in Figures [9j2 and |9j> the 
dependence of the bright tail of the pixel PDF on the in- 
strumental resolution. 9 Each curve is constructed from 100 
maps of 150 2 pixels each (after smoothing). The solid, long- 
dashed, short-dashed, and dotted curves in each panel have 
A9 = 1', 0.5', 0.24', and 0.12', respectively (or physical res- 
olutions ranging from w 110/T 1 kpc to w 13/i _1 kpc, about 
twice as large as the simulation resolution). 

8 One might naively expect a ( 1 + z)~ 3 dependence from surface bright- 
ness conservation (note that we work in units of photons cirr 2 s sr~' A -1 ). 
However, in reality the surface brightness increases slowly with redshift be- 
cause the IGM density, and hence the recombination rate, increases with red- 
shift. 

9 Decreasing A9 broadens the low surface brightness peak of the PDF but 
leaves the mean emission unaffected. 




log $ (photons cm -2 s _1 sr _1 ) 

FIG. 9. — Pixel PDFs. All curves assume z = 0.15, Az = 0.01, simu- 
lation metallicities, and the HM01 background unless otherwise specified. 
(a): The VI A1032 transition with A0 = 1', 0.5', 0.24', and 0.12' (solid, 
long-dashed, short-dashed, and dotted curves, respectively), (b): Same as 
(a), but for the C IV A1548 transition, (c): The O VI A1032 transition with 
z = 0.15, 0.10, and 0.05 (solid, dashed, and dotted curves, respectively), with 
fixed physical resolution of 25/r 1 kpc. (d): PDFs for the HMOls, HM01, 
HMOlw, and HM96 ionizing backgrounds (long-dashed, solid, short-dashed, 
and dot-dashed curves, respectively) and for no ionizing background (dotted 
curve), assuming A6 = 0.24'. 



We see that the bright tail of the PDF retains qualitatively 
the same shape, so our conclusions are essentially indepen- 
dent of angular resolution. Increased angular resolution tends 
to decrease the fraction of moderately bright pixels, by up to 
an order of magnitude, but increase the peak brightness. The 
C IV A 1548 distribution changes more with A8 than the O VI 
A1032 distribution does because C IV A1548 regions tend to 
be more compact. Resolution effects may begin to saturate at 
Af? < 0.24', but this is probably an artifact of being near to 
the simulation resolution. 

4.5. Redshift Evolution 

The redshift chosen in our Figures is relatively high; O VI 
A 1032 has nearly redshifted to 1216 A, where interference 
from Lya will become important. Equation @ suggests (and 
our numerical results confirm) that the mean surface bright- 
ness of the low-density IGM does not evolve significantly 
between z = 0-0.15. However, high surface brightness re- 
gions are not well-described by equation @ because they 
have 5 > 100. This density regime corresponds to regions 
that have broken off from the cosmic expansion; such objects 
should be treated as individual sources. Figure^ shows the 
bright end of the O VI A 1032 PDF for z = 0.05, 0. 10, and 0.15 
(dotted, dashed, and solid curves, respectively). We have fixed 
the physical resolution at each redshift to be 25/T 1 kpc, equiv- 
alent to A0 = 0.24' at z = 0.15. We see very little evolution 
with redshift: there is only a slight increase in the number of 
pixels with a fixed surface brightness as redshift decreases. 
Although an individual source has a smaller surface bright- 
ness as it moves farther from us, the star formation rate in - 
creases rapidly with redshift ( Sprinael & Hernquist 2003b), 
so there are more (and more powerful) winds at higher red- 
shifts. We find that this effect nearly balances the increasing 
cosmic distances. 
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FIG. 10.— Pixel surface brightness PDFs for z = 0.25, Az = 0.01, and 
a physical resolution of 20/T 1 kpc. The upper and lower panels show the 
PDFs for the O VI A1032 and C IV A1548 transitions. In each panel, the 
dotted, solid, and dashed curves correspond to the G4, G5, and G6 runs, 
respectively. The insets show the bright end of the PDF. All curves assume 
the HM01 ionizing background. 



4.6. Ultraviolet Background 

Figure|9j/ shows how the choice of UV background affects 
the pixel surface brightness distribution. The long-dashed, 
solid, and short-dashed curves show the PDFs for the O VI 
A 1032 transition assuming the HMOls, HM01, and HMOlw 
backgrounds, respectively. The dotted line assumes pure col- 
lisional ionization, while the dot-dashed line assumes the 
HM96 spectrum. We see that the brightest pixels are essen- 
tially independent of the choice of ionizing background. This 
is to be expected, because the peak emissivity shown in Fig- 
ure n comes from collisionally ionized gas whose location is 
nearly independent of the background choice. We note that 
the shape of the low-surface brightness peak does depend on 
the background: the distribution skews toward higher flux val- 
ues in weak ionizing backgrounds, because low-density gas is 
not photoionized to higher states (see Figure ^ an d discus- 
sion thereof). Although we do not show it, the ionizing back- 
ground also has no effect on bright C IV A 1548 emission. 

4.7. Simulation Resolution 

One important concern is how the resolution of our simu- 
lation affects the results. This is particularly important be- 
cause each bright region has only of order ~ 10 bright par- 
ticles; we therefore do not expect their structure to be fully 
resolved. Moreover, our wind prescription is stochastic and 
will not fully capture the dynamics of winds in dwarf galax- 
ies. To estimate the importance of our simulation resolution, 
we have repeated our analysis for two other simulations, G4 
and G6, which are identical to the one described above (G5) 
exc ept that they contain 144 3 and 486 3 particles, respectivel y 
(see lSpringel & Hernauistll2003bl and lNaeamine et alJl2003h . 
We show the pixel PDFs for the different simulations in Fig- 
ure m eacn panel, the G4, G5, and G6 simulations are 
shown with dotted, solid, and dashed lines. Both panels as- 
sume a uniform metallicity (fixed at that of the G5 run); we 
have therefore removed the effects of the different metal dis- 
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tribution between the simulations. We see that the simulation 
resolution has only a small effect on the resulting PDF: in par- 
ticular, the number of bright pixels increases by only a small 
amount. We thus conclude that the temperature and density 
structure of the IGM have converged to good accuracy. The 
agreement is in fact surprisingly good, given that the G6 sim- 
ulation resolves large galaxies and winds better but also adds 
many more small galaxies. 

The dependence on resolution in the simulated metallic- 
ity case is slightly larger than shown here for the bright pix- 
els, although the change at small surface brightness is non- 
negligible between the three resolutions. This suggests that 
the metal distribution does change between the simulations, 
probably because of the stochastic nature of the enrichment 
process. If each wind is composed of more particles, any in- 
dividual volume element in the simulation is more likely to 
contain such a particle, so the metal distribution is less patchy. 
We thus find that the G6 simulation has ^25% more pixels 
with $ ~ 10~ 3 -1 photons cm" 2 s" 1 sr" 1 than the G5 simula- 
tion does. The patchiness is not as significant for the bright- 
est regions (where the resolution has about the same effect 
in both the uniform and simulated metallicity cases), because 
such regions are compact and the "density" of enriched parti- 
cles is still large. Thus the stochastic nature of the enrichment 
process only has an important effect at large distances from 
the source, where the surface brightness is so low as to be 
unobservable anyway. 

4.8. Angular Correlation Functions 

We define the (flux-weighted) correlation function w{9) on 
a pixelized map as (iCroft et al.l2001l) 

wm=(5$(r)5*(r + 9)), (5) 

where r is the angular position of a pixel, (5$(r) = $(r)/$- 1, 
and $ is the mean pixel surface brightness. We measure w(9) 
for each pixel pair of the specified separation and compute the 
mean value across a given map. 

Because spectral line studies automatically contain redshift 
information, angular clustering statistics such as w(9) are per- 
haps not the optimal measure: one could use the line redshifts 
to study clustering in three dimensions. However, the angu- 
lar correlation function is considerably simpler to calculate 
for a survey geometry such as ours, and we believe it to be 
adequate for a first approach to the problem. A simple way 
to retain some redshift information is to compute w{9) only 
for sources within a relatively narrow redshift window. Shot 
noise will obviously decrease with bandwidth because more 
sources are included in the map. However, using a large band- 
width would sacrifice all redshift information. We therefore 
chose to calculate w(9) in slices like those described in the 
previous section. This preserves some redshift information 
but reduces the shot noise to a manageable level. 

Figure ITTlshows the angular correlation functions w{9) for 
the O VI A 1032 and C IV A 1548 transitions (left and right 
panels, respectively). The solid and dashed curves are the 
mean correlation function measured from fifty 2.5° x 2.5° 
slices with z = 0.15, Az = 0.01, and A9 = 1', assuming sim- 
ulation and uniform metallicities, respectively. 10 The solid 
hexagons are the median values (assuming simulation metal- 
licities) measured between maps. The error bars show the 
upper limits of the first and third quartiles; they illustrate 

10 The typical scale of emitting regions at z = 0.15 is ~ 1', so we are 
calculating the correlation function between individual systems. 
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FIG. 11.— Angular correlation function for the O VI A 1032 (left) and CIV 
A1548 (right) transitions in slices with z = 0.15, Az = 0.01, and A9 = 1'. The 
upper and lower panels show w(6) in logarithmic and linear units, respec- 
tively. The solid and dashed curves show the mean signal for the simulation 
and uniform metallicity models, respectively. The hexagonal points are the 
median values of w(9) and the error bars show the upper limits of the first 
and third quartiles (assuming simulation metallicities). The solid and open 
squares show the autocorrelation of the maps for simulation and uniform 
metallicities, respectively. 



the dispersion in measurements of w{9) between individual 
2.5° x 2.5° maps. The solid and open squares show the au- 
tocorrelation of the emission; note that the location along the 
0-axis is arbitrary because 9 = does not appear on the log- 
arithmic scale. Because the correlation function is very large 
for small separations, the top panels show w(9) on a logarith- 
mic scale. The correlation rapidly vanishes for larger separa- 
tions; in order to show this decrease, we plot w(9) on a linear 
scale in the bottom panel. 

Metal line emission is clearly clustered: w(9 = 7.5') ~ 1 for 
O VI A1032 and w(9 = 5') ~ 1 for C IV A1548, corresponding 
to physical separations of ~ 800/T 1 kpc and ~ 550/i _1 kpc, re- 
spectively. A power law of index n ~ 3.5 (n ~ 4) adequately 
fits the correlation function for O VI A1032 (C IV A1548). 
The correlation function is much steeper than that of galax- 
ies, at least partly because we have weighted the pixels by 
their fluxes. A difference from galaxies is not surprising be- 
cause of the complicated physics behind line emission: there 
is no reason to expect the line emission to correspond to bright 
galaxies. Rather, emission is likely to come from those galax- 
ies that (for whatever reason) host an ongoing wind. Note that 
the correlation length we measure is approximately the size of 
galactic winds (see Figure|4]and discussion thereof). 

The case with simulation metallicities is somewhat more 
strongly clustered at small scales because of the metallic- 
ity gradient between enriched, emitting regions and the (rel- 
atively) pristine, low-density IGM. However, differentiating 
the two through clustering will be difficult because the disper- 
sion between maps is nearly as large as the difference between 
the two models. Cosmic variance is also very large, because 
the correlation function is dominated by rare, bright pixels. 
Comparison of the mean and median shows that a typical map 
underestimates the true correlation coefficient, particularly at 
intermediate scales. This is because the positive part of the 
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signal is dominated by pairs of bright pixels, for which the 
probability is quite small. Thus, a precise measurement of 
the clustering properties of metal line emitters will require a 
survey over a large volume. 

We have not included any background emission in this cal- 
culation of w(8) because (as discussed above) the background 
level at these wavelengths is so uncertain. However, comput- 
ing the correlation function directly from the raw map (with- 
out background subtraction) smooths the correlation substan- 
tially because the total flux from a map of this size is domi- 
nated not by the metal emitters but by the background radia- 
tion field. A better strategy is to set all pixels below a certain 
signal-to-noise to zero. This procedure does not decrease the 
number of individual sources by much, because most contain 
bright cores well above the background. 

5. DISCUSSION 

We have examined the emission expected from the IGM 
in two UV metal line doublets, O VI AA1032, 1038 and C IV 
AA1548,1551. Using state-of-the-art cosmological simula- 
tions that include prescriptions for star formation and galac- 
tic winds, we computed the surface brightness distribution of 
these lines. We assume the IGM to be photoionized by a dif- 
fuse background, and we include several different models for 
this (currently uncertain) background radiation. We also con- 
sidered two different metallicity distributions for the IGM: a 
uniform metallicity model and a model with the metallicities 
predicted by our simulation (essentially including winds from 
galaxies at low and moderate redshifts). 

Can UV metal line emission from the IGM plausibly be ob- 
served? A comparison of Table 1 and Figure[8]shows that, in 
the low-density IGM, the emission from either of our tran- 
sitions is several orders of magnitude smaller than that of 
any realistic background radiation; we therefore conclude that 
these transitions provide little information about the "missing 
baryons" (i.e., baryons with high temperatures but low densi- 
ties). Instead, Figures |4] and |6] show bright emission is con- 
centrated into isolated regions with (1) significantly higher 
metallicity than the average, (2) T ~ 10 5 K (for C IV A1548) 
or T ~ 10 s 5 K (for O VI A1032), and (3) a density much 
greater than the cosmic mean. C IV A 1548 emission is gen- 
erally confined to compact regions of size < 40/T 1 kpc, while 
O VI A 1032 emission is generally more extended with bright 
central cores of size ~ 50/T 1 kpc surrounded by weaker, ir- 
regular halos. The overall space density of emitting regions is 
«ovi ~ 7 x 10~ 3 /i 3 Mpc" 3 and « C iv ~ 3 x 10~ 3 /i 3 Mpc" 3 (co- 
moving) at z = 0.15. Overall, we find that O VI A1032 has 
more promise as an observational probe. 

Detection of metal line emission is feasible (but challeng- 
ing) with existing technology. The ideal instrument for these 
observations would be an integral field spectrograph with a 
large field of view; resolution of at least ~ 10" would be nec- 
essary to study the detailed structure of the sources. Existing 
instruments do not meet these criteria, though the next gener- 
ation of instruments probably will. The Galaxy Evolution Ex- 
plorer (GALEX), with a large field of view but relatively low 
spectral resolution (~ 10 A), may detect the brightest cores, 
though its poor spectral resolution would make it impossi- 
ble to identify doublets unambiguously. A project underway 
to construct a balloon-borne wide-field UV spectrograph will 
have a limiting sensitivity $ ~ 500 photons cm" 2 s" 1 sr" 1 in a 
single night's observation (D. Schiminovich, private commu- 
nication); such an instrument could detect the cores of many 
O VI emitting regions. Finally, the proposed Spectroscopy and 



Photometry of the IGM's Diffuse Radiation 11 would have the 
sensitivity to detect the bright cores. Background-limited ob- 
servations to detect emission on ~ 100 kpc scales will require 
a much larger collecting area than currently planned for any 
wide-field instrument. High spectral resolution would help 
to reduce the background, which is especially important for 
observations redward of 1216 A. 

Comparison of the topology of luminous regions to the 
temperature and metallicity distribution in the IGM (see Fig- 
ure |4j shows that most bright particles are embedded in hot 
(T > 10 6 K) enriched regions characteristic of galaxy groups 
and winds. We therefore suggest that these transitions probe 
models of metal enrichment and galactic feedback. We have 
shown that these observations can distinguish a uniform IGM 
metallicity from one in which the metals are dispersed by 
powerful winds. If mixing is more efficient than the simu- 
lation allows, the emission could be more extended than we 
predict but have smaller peak luminosities. Unfortunately, the 
feedback parameters in our simulation model are fixed, and 
we are thus unable to investigate how the emission varies with 
them. We have deliberately chosen an extreme wind model (in 
which none of the supernova energy is lost to radiative cool- 
ing before powering the wind). Metal emission will tend to 
be more compact if winds have less energy simply because 
they cannot penetrate as much of the IGM. However, a more 
subtle consequence of our choice is that it allows all of the 
supernova energy to heat the wind medium. Weaker winds 
may therefore lead to more emission, because the wind stays 
cooler near to the host galaxy; however, note that th e wind 
veloc ity of our model is consistent with observations jMartinl 

1999) . 

It is important to note that the extent of observable O VI 
A 1032 and (especially) C IV A 1548 emission does not trace 
out the extent of metal enrichment. Collisional excitation of 
these lines requires high density and T > 10 5 K. As winds 
expand into the IGM, their density decreases rapidly and ther- 
mal energy injected by the wind is lost to cooling. Thus we 
would expect emission from wind-enriched regions to fade 
with time. The declining star formation rate at low redshifts 
implies that active winds become relatively rare in the local 
universe. Most enrichment occurred at earlier epochs, giv- 
ing the winds time to expand and mix with the surrounding 
IGM. Nevertheless, observations of line emission will provide 
a census of enriched, warm regions at the present day. 

Higher resolution observations of individual sources have 
the potential to teach us about the structure of galactic winds. 
If viewed in O VI A 1032 emission, we predict a central, 
bright core surrounded by an irregular, low surface brightness 
halo. However, the actual structure will be much more com- 
plex than our simulation, with its relatively coarse resolution, 
can accurately describe. Real galactic winds have extremely 
complicated internal structures with several different temper- 
ature phases that can be separated through observations of 
diffus e X-ray gas with kT > 0.1 keV (e.g., IStrickland et alJ 

2000) , absorption lines of low-ionization species (such as 
Na I;|H eckman et al. 2000), and bot h emission lines (like Ha; 
Martinl [l998l) and absorption lines ( Shaplev et al. 120031) typ- 
ical of highly ionized regions. These are thought to probe, 
respectively, a hot phase with large filling factor, a cool phase 
made up of embedded cl ouds, and a warm phase on the in- 
terface between the two ( Heckmanl 120021) . This last phase 
should host O VI A 1032 and C IV A 1548 emission. Re- 

1 1 See http://www.bu.edu/spidr/indextoo.html 
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cent observations of the wind surrounding M82, a nearby 
starburst galaxy, show that O VI A 1032 is not an efficient 
coolan t of gas withi n several kiloparsecs of the host galaxy 
jHoooes et alJ 120031) : however, the limiting surface bright- 
ness (<J> < 10 5 photons cm -2 s" 1 sr" 1 for this nearby system) 
is compatible with our predictions. Because of the many 
observational mysteries, metal lines offer a powerful probe 
of galactic winds complementary to existing techniques. In 
particular, emission line observations, which can in princi- 
ple probe scales up to ~ 100 kpc, will help to determine how 
winds evolve on scales large compared to their host galax- 
ies. Correlating the line emission maps to galaxy surveys will 
teach us about the properties of galaxies that are able to drive 
powerful winds and by extension the types of galaxies respon- 
sible for enriching the IGM with metals. 

Our predictions are subject to several important uncer- 
tainties. First, we have excluded self-shielded and star- 
forming gas from our analysis. Self-shielded gas is un- 
likely to contribute significantly to O VI A 1032 emission, be- 
cause it is much too cool, but it could add a small number 
of additional C IV A1548 sources. Self-shielding is signifi- 
cantly more importan t for studies of hydrogen line transitions 
(Furlanetto et al.l l2003l) . because dense, cool clumps of gas 
show strong Lya and Ly/3 lines. A high enough spectral reso- 
lution to resolve the O VI and C IV doublets will be necessary 
to prevent contamination. Galaxies will inevitably be strong 
emitters at these wavelengths and present another potential 
source of confusion for studies of the IGM. However, they 
can be removed from maps like the ones we show because 
they will also be strong continuum emitters, provided that the 
angular resolution of the instrument is smaller than the size of 
the emitting regions. 

Another set of caveats has to do with the simulation itself. 
The mo st obvious is its finite resolution. However, we showed 
in 34.7l that the density and temperature structure of our sim- 
ulation has converged. The remaining uncertainty is in the 
patchiness of the metal distribution, but that has little effect 
on the bright regions. Another concern is that the simulation 
neglected metal line cooling. We argued in 33. II that this is 
probably only important in highly enriched regions such as 
winds. We expect it to decrease the number density of bright 
regions by, at most, a factor of a few, and to increase their sur- 
face brightness (because they cool thro ugh the relevan t tem- 
perature regime sooner). Note that Hoopes et al. (2003) argue 
that line cooling is not important for winds within ~ 10 kpc 



of the host galaxy. 

Finally, we neglected local effects such as dust and vari- 
ations in the background radiation field. Because luminous 
regions tend to be quite dense, the local radiation field likely 
will not strongly affect our results. Most important, we find 
that the bright end of the pixel distribution is largely unaf- 
fected by uncertainty in the ionizing background (see Fig- 
ure |9jf), although ionizing radiation from the host galaxy of 
a wind could still have substantial effects. Dust is a particular 
concern because galactic winds have been ob served to carry 
it over sizable distances (Heckman et al. 2000); much of our 
emission comes from regions recently enriched by winds, so 
dust will likely decrease the emissivity of such regions. For- 
tunately the measured extinction in local starburst winds is 
usually not large and may be patchy, so we expect our quali- 
tative results to remain unchanged. 

Despite these uncertainties, our results show that C IV 
A 1548 and especially O VI A 1032 emission can be detected 
from the IGM. While emission from the low density "missing 
baryons" is too weak to be detected, we find that observations 
of these transitions may distinguish different enrichment pat- 
terns in the IGM; they thus offer one of the few methods to 
probe the effects of mechanical, wind-driven feedback from 
galaxies. Many of the uncertainties we have mentioned have 
to do with the detailed structure and mechanism of winds; ob- 
servations of the metal emission will of course help to resolve 
these questions. When combined with observations of emis- 
sion from higher ion ization states (such as O VII and O VIII; 
Yoshikawa et al J 120031) . maps of the UV doublets we have 
studied will also teach us about the thermal structure of the 
enriched IGM, offering an even more complete picture of the 
wind mechanism. 
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